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Abstract 

We investigate the A{tt^ , K^)j>^A reaction within a covariant modeL We con- 
sider those amphtudes which are described by creation, propagation and decay 
into relevant channel of iV*(1650), 7V*(1710), and iV*(1720) intermediate bary- 
onic resonance states in the initial collision of the incoming pion with one of 
the target nucleons. The bound state nucleon and hyperon wave functions are 
obtained by solving the Dirac equation with appropriate scalar and vector po- 
tentials. Expressions for the reaction amplitudes are derived taking continuum 
particle wave function in the plane wave approximation. Numerical calculations 
are presented for reactions on ^■^C, ''"Ca, ^^V and ^^Y target nuclei. The predic- 
tions of our model are in reasonable agreement with the available experimental 
data. 

Keywords: strangeness production, pion- nucleus collisions, covariant model 
PACS: 25.40.Ve, 13.75.-n, 13.75.Jz 



1. Introduction 

Hypernuclei, where one or two nucleons (N) are replaced by hyperons {Y) 
in the bound orbits, provide a unique opportunity for studying a new form 
of the hadronic system which has strangeness degrees of freedom [E 0, 0, 01 ■ 
Lambda (A) hypernuclei are the most familiar and extensively investigated hy- 
pernuclear systems. Since the A hyperon does not suffer from Pauli blocking 
by other nucleons, it can penetrate deep inside the nucleus and form deeply 
bound hypernuclear states. Thus, hypernuclei can provide information about 
the nuclear states which are not accessible in ordinary nuclei. Such systems are 
perhaps the only tool currently available to get information about the A- TV in- 
teraction as A-nucleon scattering experiments are very difficult to perform due 
to short life time of the A particle. During the past years, data on hypernu- 
clear spectroscopy have been used extensively to extract information about the 
hyperon-nucleon interaction within a variety of theoretical approaches (see, e.g., 

Refs. [ilSSBiEliil,!!!). 



A hypernuclei have been studied extensively by the stopped as well as the 
in-flight {K~,TT^) reaction (see, e.g., the reviews @, 0, 0, 13, 14, 15|) and also 



(see, e.g 


16|, 


17 





by the (7r+,if+) reaction [16|, [iTl, [iSl, ll9| . The kinematical properties of the 
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{K~,Tr~) reaction allow only a small momentum transfer to the nucleus (at 
forward angles), thus there is a large probability of populating A-substitutional 
states (A assumes the same orbital angular momentum as that of the neutron 
being replaced by it). On the other hand, in the (7r+,iir+) reaction the mo- 
mentum transfer is larger than the nuclear Fermi momentum. Therefore, this 
reaction can populate states with the configuration of an outer neutron hole and 
a A hyperon in a series of orbits covering all the bound states having high spin 
natural parity configurations. The richness of the spectroscopic information on 
A bound states in the (tt^, K^) reaction was demonstrated in the experiments 
performed at the Brookhaven National Lab and National Laboratory for High 
Energy Physics (KEK) (see, e.g., Ref. [1] for a comprehensive review). Fur- 
thermore, although the reaction cross section of the strangeness production, via 
the {tt^,K^) process, is smaller than that of the strangeness exchange reac- 
tion (K^,TT^), the higher luminosity of pion beams makes experiments more 
feasible. 

In the experimental studies reported in Refs. [H, this reaction has been 
used to carry out the spectroscopic investigations of hypernuclei ranging from 
light mass ]^C, to medium mass ^^V and ^^Y with the best resolution (^1.6- 
1.7 MeV) achieved in the spectrometer at KEK. This experiment has succeeded 
in clearly observing a characteristic fine structure in heavy systems by precisely 
obtaining a series of A single-particle states in a wide range of excitation energies. 

Most of the theoretical models used so far to describe the (7r+,if+) reac- 
tion employ a non-relativistic distorted wave impulse approximation (DWIA) 
framework [i^l (see also Ref. Q for a comprehensive review of these models). In 
these calculations, the A bound states are generated by solving the Schrodinger 
equation with Woods-Saxon or harmonic oscillator potentials. However, for 
processes involving momentum transfers of typically 300MeV/c or more, a non- 
relativistic treatment of the corresponding wave functions may not be adequate 
as in this region the lower component of the Dirac spinor is no longer negligible 
in comparison to its upper component (see, e.g., Ref. (2ll|). 

In this paper, we study the A(7r"'", K'^)/^A reaction within a fully covariant 
model by retaining the field theoretical structure of the interaction vertices and 
by treating the baryons as Dirac particles. In this model, the kaon production 
proceeds via the collision of the projectile pion with one of the target nucle- 
ons. This excites intermediate baryon resonance states (N*) which decay into 
a kaon and a A hyperon. The hyperon is captured in the respective nuclear 
orbit while the kaon rescatters onto its mass shell (see Fig. [T]) . A similar pic- 
ture has been used to describe the A{p, K'^)\B and A{'j, K^)\B reactions in 
Refs. [l^, [2^, [131 . In our model, the intermediate resonance states included are 
N*{1650)[^~], iV*(1710)[i^], and A^* (1 720) [f_^ which have dominant branch- 
ing ratios for the decay to the fT+A channel [1^, [1^. Terms corresponding to 
the interference among various resonance excitations are included in the total 
reaction amplitude. 

In section [21 we present the details of our formalism for calculating the 
amplitudes corresponding to the diagrams shown in figure [H In section [SI 
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Figure 1: Graphical representation of our model to describe the (7r+,_fS'+) reaction. The 
elliptic shaded area represent the optical model interactions in the incoming and outgoing 
channels. 



numerical results are presented for the (7r+, K^) reaction on ^^C, ""^Ca, ^^V and 
^^Y targets using continuum wave functions in the plane wave approximation. 
Summary, conclusions and future outlook of our work are given in section [6l 

2. Covariant model for the A{'7t~^ , K~^)aA reaction 

The structure of our model for the (7r+, K^) reaction is described in Fig. [TJ 
The A^* corresponds to the N*{m50)[^^], iV*(1710)[i^], and TV* (1720) [f^] 
baryon resonance intermediate states. Terms corresponding to interference be- 
tween various amplitudes are retained. The elementary process involved in this 
reaction is shown in Fig. [2 

It is clear that our model has only s-channel resonance contributions. In prin- 
ciple, Born terms and the resonance contributions in u- and t-channels should 
also be included in description of both the processes depicted by Figs. [1] and [H 
These graphs constitute the non-resonant background terms. Their magnitudes 
depend on particular models used to calculate them and also the parameters 
used in those models. This can be seen in Ref. ^] where the effect of the 
background terms is studied for 5*11 phase shifts in the pion-nuclcon interac- 
tion. Whereas the contributions of the background terms are about 15-20 % 
of the resonance terms within a coupled-channel K matrix model, they are 



limited to less than 10% in the model of Ref. [28[ for the energies of our inter- 
est (« 1 GeV/nucleon; the corresponding invariant mass is about 1.7 GeV) (see 
Fig. 15 of [27|). The effect of the background terms on the total production 
cross sections of the ir^p — > K^A reaction at our beam energies can be indi- 
rectly inferred from the calculations reported in Ref. ^2^ . It is seen in Fig. 2 of 
this reference that the total production cross sections at invariant mass around 
1.7 GcV are dominated by the contributions of 5*11 and Pn resonance terms. 
Therefore, the magnitudes of the background terms are likely to be limited to 
about 10-20 % of those of the resonance terms. Our results reported later on 
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Figure 2: Tree diagram for the elementary process of pion-induced strangeness production via 
resonance excitation and decay on a single nucleon. 



in this paper may be uncertain by 10-20 % due to omission of the background 
contributions. 

2.1. Interaction Lagrangians 

For the interaction terms of the spin-1/2 resonances, we have vertices of 
pseudoscalar (PS) or pseudovector (PV) form. The pseudovector couphng is 
consistent with the chiral symmetry requirement of the fundamental theory of 
strong interactions (quantum chromodynamics (QCD)). In contrast to that, the 
pseudoscalar one does not have this property, but it is easier to calculate. The 
couplings are in both cases fixed in such a way that they are equal on-shell; 
for off-shell cases, their difference is suppressed due to the denominator of the 
resonance propagator. It is, therefore, arguable which Lagrangian to use. The 
best approach would be to introduce a mixing parameter, which was investigated 
m Refs. M,\ 



To avoid the introduction of additional parameters in our model due to 
a PS-PV mixing for the interaction Lagrangians of Refs. 3^, 24|, we use the 



convention of either choosing the PS or the PV couplings for these vertices as 
done in Refs. 27|, |31|. This is in line with the studies reported in Refs. [13, n 



The pseudoscalar interaction Lagrangians for the spin-1/2 resonances are 
given by 

^Inn' . = -9iiNN'ipN''r{'r ■ 4>^)'>Ijn + h. c. , (la) 



1/2 



^N', KA = -gN'AKlpN'TcjjKi^A -|- h. C. , (lb) 

1/2 

where the F takes care of parity conservation. We use 

1 for odd parity 



F = 



ij^ for even parity 



and h. c. in Eqs. ([T]) denotes the hermitian conjugate. 

The pseudovector Lagrangians involve the derivative of the pion wave func- 
tion rather than the wave function itself. This introduces an additional mass 
dimension, which is taken care of by a "rescaling" of the coupling constant. It 
also ensures the matching of the on-shell behaviour the two types of Lagrangians. 
The pseudovector Lagrangians are given by 

= ^-^^^i^N'rrd^{T-4>Ji^N + h.c., (2a) 

£PY ^ 9N'KA ^^,^Mra^^^^^ + c. , (2b) 

1/2 mjv* ± ITT-A 
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Table 1: Coupling constants for various vertices used in the calculations. 



vertex 


coupling constant (g) 


iV*(1650)A^7r 


0.81 


N*{1650)AK+ 


0.76 


iV*(1710)A^7r 


1.04 


N*{1710)AK+ 


6.12 


N*{1720)Nn 


0.21 


N*{1720)AK+ 


0.87 



where T is given by 

r 



i for odd parity 
7^ for even parity 



and the upper and lower signs are used for even and odd parity resonances, 
respectively. 

The spin-3/2 resonance Lagrangians are given by 

^irWAf- , = '^""^^ ■ip%.df,{T ■ (p^jtpN + h. c. , (3a) 

^n;,^ka = '^^ ^^ 'ipN'^^fj.'I^K^A + h- c. . (3b) 

Such a form of the coup ling was used for studying the hypernuclear production 
also in other reactions |2ll. [25t. The values and signs of the various coupling 
constants have been taken from Refs. and are shown in table [1] These 



parameters describe well the associated K'^A production in proton-proton col- 
lisions within a similar resonance picture. All the pion-resonance-kaon vertices 
that are of interest in this paper are involved in this reaction. Thus the ver- 
tex parameters used by us inherently describe the elementary process shown in 
Fig.H 

2.2. Resonance propagators 

The two interaction vertices of figure [5] are connected by a resonance prop- 
agator. For the spin- 1/2 and spin-3/2 resonances the propagators are given 

by 

_ j^p" + m 

- V-(m-ir^./2)2 (4) 

and 



^3/2- y _(m-irw./2)2^ ' 

respectively. In Eq. [5] we have defined 



3 3m^ 6ra 
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rjv* in Eq. (g]) and Eq. ^ is the total width of the resonance. It is introduced 
in the denominator term to account for the finite Hfe time of the resonances for 
decays into various channels. This is a function of the centre of mass momentum 
of the decay channel, and it is taken to be the sum of the widths for pion and 
rho decay (the other decay channels are considered only implicitly by adding 
their branching ratios to that of the pion channel) [2^ . We have not introduced 
any correction in the resonance propagators to account for the nuclear medium 
effects as no major chan ge i s expected in our results due to these effects. As 



pointed out in Refs. [3J, |35|, the medium correction effects on the widths of 



the s- and p-wave resonances, which make the dominant contribution to the 
cross sections investigated here, are not substantial. The reason for this is that 
resonances occur only as intermediate states which implies an integration over 
their respective spectral distributions. 



3. Nuclear model 

The spinors for the final bound hypernuclear state (corresponding to momen- 
tum Pa) and for the intermediate nucleonic state (corresponding to momenta 
Pn) are required to perform numerical calculations of various amplitudes. We 
assume these states to be of pure-single particle or single-hole configurations 
with the core remaining inert. In experimental measurements, however, core 
excited states have also been detected (see, e.g., [4j). A covariant description 
of the core polarisation can, in principle, be achieved by following the method 
discussed, e.g., in Ref. [1^1 • This procedure is somewhat tedious and is out of 
the scope of our present study. Therefore, in this paper we concentrate on those 
transitions which involve pure single-particle and single-hole states. 

The spinors in momentum space are obtained by the Fourier transformation 
of the corresponding coordinate space spinors which are solutions of the Dirac 
equation with potential fields consisting of an attractive scalar part (Vg) and 
a repulsive vector part (Vv) having a Woods-Saxon form. This choice appears 



justified as the Dirac-Hartree-Fock calculations in Refs. [37|, |38| suggest that 
these potentials tend to follow the nuclear shape. The same potential form has 
also been used in the relativistic one-nucleon model [3^ |40| | and two-nucleon 
model calculations [2l[ of the (p, tt) reaction. 

3.1. Nucleon bound states 

In our approach for describing the nucleon and hyperon wave functions, we 
begin with the Dirac equation which is modified by introducing a scalar and 
a time-like vector potential. The modified Dirac equation with potentials is 
written as 

(i7^9^ - m - - V;) ^ = . (7) 
The solution, transformed to momentum space, can be written as in Ref. [23| 



m = S{po - E) ^ ' (8) 
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Table 2: Potential parameters of the nuclear vector and scalar potentials. For C and "Ca, 
the potentials are fitted to the experimental data on the charge radius, the nucleon separation 
energy, and the first diffraction minimum of the charge form factor for the nucleon (Ref. [4l|). 
For ^'^V and *^Y, the potential depths are fitted to the neutron separation energies with radial 
parameters being the same as those of ^''Ca. 



nucleus 


K [MeV] 


f Ov [fm] 


Cv [fm] 


V; [MeV] 


ros [fm] 


Us [fm] 


12c 


385.7 


1.056 


0.427 


-470.4 


1.056 


0.447 




348.1 


1.149 


0.476 


-424.5 


1.149 


0.506 


Sly 


309.7 


1.149 


0.476 


-382.3 


1.149 


0.506 


89y 


317.8 


1.149 


0.476 


-392.3 


1.149 


0.506 



where 'p' is the four momentum vector of the particle, 'p' represents the corre- 
sponding three momentum; its magnitude \p\ is denoted by k and its direction 
by p. Since we work in momentum space throughout this work, we will drop 
the hats from the Fourier-transformed functions ^, /, and g. 
The spin-spherical harmonics yj^^ are given by 

yeT'i^)^ ^ {^,TTT'i:S,ms\j,mj)YimAx)Xs,m, , (9) 

where x is the usual two-dimensional Pauli-spinor and Yim are the spherical 
harmonics of the first kind, £' = 2j — i, and me + iris = ^j, and both are 
coupled to good total angular momentum. 

The potentials in Eq. ([7]) have radial shapes of the Woods-Saxon type. The 
depths of the potential Vba, the radius Ra = roaA^/"^, and the diffuseness Oq are 
treated as free parameters which are fixed by fitting to the experimental data 
of the charge radius, the nucleon separation energy, and the first diffraction 



minimum of the charge form factor for the nucleon (see Ref. [4l|). Table [2] lists 
the potential parameters used for the nuclear bound states. 

3.2. Hyperon bound states 

The A, in contrast to a nucleon, can occupy any bound state since it is not 
subject to the Pauli exclusion principle with respect to the nucleons. This prop- 
erty makes it an excellent probe for single bound states in a nuclear potential. 

In this case too, the spinors in the momentum space were obtained by Fourier 
transforming the corresponding coordinate space spinors which are solutions of 
the Dirac equation with potential fields consisting of an attractive scalar part 
and a repulsive vector part having a Woods-Saxon form. With a fixed set of the 
geometry parameters (reduced radii rg and and diffusenesses Og and a^), the 
depths of the potentials were searched in order to reproduce the A separation en- 
ergy of the particular state (the corresponding values are given in table[3]for 
and ^\V, and in tabled] for ^%Y). We use the same geometry for the scalar and 
vector potentials. The depths of potentials Vs and Vv were further constrained 
by requiring that their ratios are equal to —0.81 as suggested in Ref. [i^l- Ex- 
perimental inputs have been used for the configurations and the single baryon 
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Figure 3: The ^^C bound states in momentum space. 



Table 3: Potential parameters of the A vector and scalar potentials of the ^^C and the 
hypernuclei. The ^^C geometry is fixed to rov = 1.1486 fm, aov = 0.3960 fm and ros = 
1.1207 fm, aos = 0.4764 fm. The geometry is fixed to rov = 0.9827 fm, aov = 0.5779 fm 
and ros = 0.9825 fm, gps = 0.6064 fm. 



orbital 


^bind [MeV] 


K,v [MeV] 


1^0. [MeV] 


'iC(p3/2) 


10.79 ±0.11 
0.10 ±0.04 


171.5230 
171.5230 


-211.7654 
-211.7654 


'lV(si/2) 
'AV(p3/2) 
'AV(d5/2) 


19.75 
11.75 
3.75 


151.2006 
171.2202 
197.5942 


-186.6674 
-211.3830 
-243.9435 



binding energies of the states in hypernuclei ^aC in Refs. [3, [l^, 14 1, ^aV in 
Ref. 3, and ^^Y in Ref. Q. However, the single baryon binding energies for 
states in *XCa hypcrnucleus were taken from the density dependent relativistic 
hadron field (DDRH) theory predictions of Refs. [HI, [l3li which reproduce the 
corresponding experimental binding energies reasonably well. In this case, we 
have compared, for each state, the spinors calculated by our well-depth search 
method with those calculated within the DDRH theory (see Refs. [Ill, [l2|) and 
find an excellent agreement between the two. 

The upper and lower components of the Dirac spinors for the A bound states 
in momentum space are shown in figures [SHS] for ^^C, ^AV, and *^Y, respectively. 
We note that in each case, only for momenta q < 2 fm^^ is the lower component 
of the spinor substantially smaller than the upper component. In the region 
of momentum transfer pertinent to exclusive kaon production in pion-nucleus 
collisions, the lower components of the spinors are not negligible as compared 
to the upper component which clearly demonstrates that a fully relativistic 
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Table 4: Potential parameters of the A vector and scalar potentials of the hypernucleus. 
The geometry is fixed to roy = 0.9827 fm. (uk = 0,o779fm and ro^ = 0.9825 fm, aos = 
0.6064 fm. 



orbital 


Ebind [MeV] 


K)v [MeV] 


Vos [MeV] 




23.1 


149.3054 


-184.3276 




16.5 


171.1863 


-211.3411 


'RY(p3/2) 


16.5 


165.3949 


-204.1913 




10.0 


208.0682 


-256.8744 




10.0 


190.7247 


-235.4626 


'iY(/5/2) 


2.3 


253.8945 


-313.454 




2.3 


214.0216 


-264.2242 



CO r, 

s;- 10-1 

CO n 

^ w--> 

^ 10-' 



12 3 4 

q [fm-i] 

Figure 4: The ®XV bound states in momentum space. 




approach is better for an accurate description of this reaction. 



4. Calculation of the cross section 

Wc use the subscripts n, K, A, and B to denote the quantities of the in- 
coming pion, the outgoing kaon, the target nucleus, and the final hypernucleus, 
respectively. The differential cross section for the (tt"*", K~^) reaction is given by 



da = 



1 d pk d pb mATUB 



(27r)2 2Ek 2Eb \pJV^ 



S^^HPn+PA-{PK+PB)) . (10) 



The summation is carried out over initial (mj) and final (m/) spin states; J2r 
indicates the summation over all three resonances. 
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Figure 5: The bound states in momentum space. 



After having established the effective Lagrangians and the couphng con- 
stants, one can write down, by following the well known Feynman rules, the 
amplitudes for the graph shown in figure O The isospin part is treated sepa- 
rately which gives rise to a constant factor for each graph. Thus, the matrix 
element for our process is given by 



M 



d^fcAT 



d-^/cA f d'^p 



(27r)4 J (27r)4 J (27r)4 



X 1 



TUN' 



:Ti3(t>T,{p - kN)il'N{kN) ■ (11) 



In Eq. pip , the factors and F^ are given by the interaction Lagrangians 
from Eqs. J!])-©. The (/)s denotes the meson wave functions, and the -08 are the 
solutions of the in- medium single-particle Dirac equations given in sections 13.11 
andO 

The incident pion and outgoing kaon fields are given by 



{-lY''Ye^m^{pK)Y;^^Jp'K) 



P^k^*{p'k)^S{p'j,, 



E 



Ek) 

X fiK{k'K-,kK) , 



(12) 
(13) 
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where i?^ and Ek represents the energies of the incident pion and outgoing kaon, 
respectively; px and p^^ denote the meson on-shell momenta. The functions fi 
are given by 

oc 

Mk') = ^ / Ji{k'r)ffirydr , (14) 



where the wave function is the coordinate space sohition of the Klein-Gordon 
equation with a meson- nucleus optical potential (see, e.g., Refs. H, 3]). 



It should be mentioned here that due to the oscillatory nature of the wave 
functions in the asymptotic region, the integrals involved in Eq. (fT^ and in 
Eq. ()13p , and in a similar fashion the radial components of the bound state wave 
functions in Eq. ([H|), converge poorly. Such integrals can, however, be calculated 
very accurately by using a contour integration method as in Ref. [Zst . 

In the plane wave approximation, the wave functions (/jtt and 0]^, are 
given by 

</)W(p;) = 54(p;-p.), (15) 

<t>^K^*ip'K)^SHp'K~PK). (16) 

Thus, the integration over fcjv and k\ become redundant. This not only reduces 
the dimensionality of the integrations by a factor of eight, but also removes the 
requirement of partial wave summations altogether. 

For a fully dynamical description of the production of hypernuclei, however, 
we have to consider the interactions of the incoming and outgoing particles with 
the nucleus. In this exploratory study, however, we would like to make calcu- 
lations for very many cases involving a variety of targets and beam energies in 
order to understand the basic mechanism of this reaction and to get the relative 
estimates of the order of magnitudes of various cross sections within a covari- 
ant model. Therefore, we restrict ourselves to the plane wave treatment of the 
scattering states in the initial and final channels. The weak mutual interaction 
(kaon-nucleus) in the final channel and higher energies of the projectile in the 
initial channels do provide support to this choice. Indeed, it has been shown 
in Ref. [46?| that in the non-relativistic impulse approximation calculations the 
pattern of the spectra look similar in the PW and DW calculations — only the 
absolute magnitudes of the cross sections are affected by the initial and final 
channel meson-nucleus interactions. Indeed, distortion effects from an optical 
potential [/opt should contribute in the leading order by 0{\Uopt/ E\), which is 
well below unity at the energies around 1.0 GeV considered in this work. 

Nevertheless, we have made an estimate of the influence of the initial state 
interaction (pion- nucleus channel) on the magnitudes of the {tt~^,K~^) cross 
sections within an eikonal approximation. 

4-.1. The eikonal approximation for the pion-nucleus interaction 



The eikonal approximation (see Refs. [47|, |48[) has been quite successful in 



describing the scattering of pions on nuclei at higher incident energies [49l [& 
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In this approximation, the relative motion wave function is given by (see for 
example Ref. (isll) 



(x) = exp < ifc • a; — - J V{b, z') dz' > , 



(17) 



where k is the incident (asymptotic) momentum of the particle, v — \k\/m is 
the magnitude of the incident velocity, and b is the (two-dimensional) impact 
vector in cylindrical coordinates. The optical potential is, in general, a complex 
function V ~ U — iW which results in a phase factor from the real part U and 
an amplitude reduction from the imaginary part W. 

Since at higher energies several effects are suppressed, a simple form for 
the optical potential can be used. For example, the tgi-approximation (see, 
e.g., Ref. [48|) relates the potential to the free-space single-particle scattering 
amplitudes (or to the total cross section) and the density. The optical potential 
in this case is given by 



V;Dpt(a;) 



47r 

2i?lab 
k 

2i?lab 



[fmpQpi,'^^ "t" /'mnPn('^)] 



g{x) , (18) 



where fmp and fmn are the elementary free-space meson-proton and meson- 
neutron scattering amplitudes, respectively. Applying the optical theorem, they 
can be substituted by the total cross sections crj^* and cr^^. The respective 
ratios of the real to the imaginary part of the scattering amplitudes are denoted 
by 7mp := '^fmp/^fmp and 7,„„ := '^fmn/^fmn- We have also separated the 
neutron and proton contributions to the potential, whereby Z is the proton 
number, N is the neutron number, and A is the total number of nucleons. 

For spherically symmetric nuclei the density, and hence the potential, de- 
pends only on the magnitude of a;, r := |a;|. In these cases we can integrate 
Eq. (|17p by rewriting the argument of the potential as 



V{r) = V{b, z)^V (v4^FT^) 



(19) 



We approximate the density by relatively simple parameterisations which are 
easier and faster to compute numerically. For light nuclei [A < 16), a modified 
Gaussian shape has been used, 



1 



0G(r) 



2{A-A) r2 
3 i?^ 



(20) 



with Rg being the radial parameter. For heavier nuclei {A > 16), a Woods- 
Saxon shaped density with the radial parameter R and the diffuseness parameter 
a has been used, 

£•0 



gwsir) 



l + exp{r^} 



(21) 
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Table 5: Density parameters for the eikonal approximation for ^^C and *XCa, fitted to the 
nucleon wave functions and the elastic scattering cross sections. 





A 


Z 


Woods-Saxon 


Gauss 










i?[fm] 


a [fm] 


Rg [fm] 




12 


6 


2.24 


0.46 


1.6 




40 


20 


3.49 


0.547 


2.08 



The densities are normalised to the total nucleon number, such that 

j q{x) d^x^A, (22) 

which is already fulfilled in the Gaussian case. In the Woods-Saxon case this 



determines qq, which is given by Refs. [5lL |5 



3^ 1 

^°^4^i?3i + (-)2 ■ (23) 

We fit the parameters Rq, R and a to the radial densities. The fitted values 
for both density approximations for ^^C and ^XCa are given in table El 



5. Results and discussion 

5.1. Cross sections for the A{tt^ ,K^){^A reaction 

First, we present our results for the ^^C(7r+, K^Y\G reaction at an incident 
pion lab momentum (piab) of 1.05 GeV. In the left panel of Fig. [HI we show a 
comparison of our calculations with the experimental data (taken from Ref. [1^ ) 
for the angular distribution of the ^^C(7r+, K^Y\G reaction where the hy- 
pernuclear state has the [(P;^/2)^' ('^i/2)a] configuration. The leading contribu- 
tion to this transition comes from the 1" ground state of ^aC corresponding 
to the first peak seen in the hypernucleus spectrum. Also shown in this figure 
are the contributions of the individual baryon resonances. We note that while 
the contributions of the TV* (1650) and iV*(1710) resonances are almost iden- 
tical, those of the TV* (1720) state is weaker by factors of 3-5. The dominant 
contribution of the two spin- 1/2 resonances was also noted in case of the hy- 
pernuclear production reactions studied via {■k~^,K'^) and {'^,K~^) reactions in 



Refs. [22, |53] within a model similar to that employed in this paper. We also 
note that interference effects of the resonances are important as their individual 
contributions do not sum up to the total cross sections shown by the solid lines. 
We see that the agreement between the calculations and the experimental data 
is better at forward angles while at higher angles the calculation overestimate 
the data. This trend was also noted in the non-relativistic distorted wave im- 
pulse approximation (DWIA) calculations as shown in Ref. [l^. This indicates 
that at larger momentum transfers the simple no-core excitation picture may 
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Figure 6: The differential cross section for the 7r++^^C — > if + +^^C reaction leading to states 
of the hypernucleus with [{p^/\)n , {1^1/2) a] (left panel) and [(P3/2)iV> {P3/2)a] (right 

panel) configurations. The lab momentum of the incoming pion (piab) is 1.05 GeV. 



not be adequate. Indeed, in Ref. [3] it is shown that agreement of the DWIA 
calculations with the data improves if contributions from the core excited states 
are also included in the calculated cross sections. There could also be a weak 
contribution (about 10%) from the 2~ member of the multiplet. 

The next A-bound state in ^aC has the configuration [(^"^^2)^7 (p)a] which 
according to Ref. [l3| is just bound by about 0.1 MeV. In the right panel of Fig. El 
we compare the calculated differential cross section for this transition with the 
data (set 5 of Ref. [loj). The theoretical results correspond to the natural parity 
2+ state of the hypernucleus. In this case our calculations underestimate the 
data at very forward angles. It should, however, be noted that the data could 
also have contributions from other members of the [(^1^/2)-'^' (p)a-] configuration 
and also from the core excited configurations. Nevertheless, the DWIA calcu- 
lations, which include contributions from such configurations, overestimate the 
data for this transition 19] in the entire angular range. 

In Fig. [7] we show the contributions of various baryonic resonances to the 
angular distribution of the 7r++*°Ca — > K^+^%Ca. reaction corresponding to the 
2+ state of ^ACa having a [(^3/2)^1 (si/2)a] configuration. The corresponding 
binding energy was taken to be 20 MeV which is consistent with the experimental 



value reported in Ref. 5j|. We note that like the ^^C, case the two spin-1/2 



resonances contribute almost equally to the cross sections and both are larger 
than the contribution of the spin-3/2 one. Detailed experimental data are not 
available for this case. The single point reported in Ref. [H3| corresponds to the 
differential cross section at zero angle which is ^9mb/sr. Our result for this 
case is close to this value. 

In Fig. [5] we show the differential cross section for the same reaction as that 
in the left panel of Fig.Elas function of the magnitude of the momentum transfer 
to the nucleus {q — — pj^) for several beam energies. We see that shapes 
of the differential cross sections are qualitatively different from each other at 
different beam energies. However, their absolute magnitudes differ at forward 
angles which is more prominent as the beam energy is increased from 1.0 GeV 
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to 1.4 GeV. However, for beam energies > 1.4 GeV this difference is small. At 
larger momentum transfers the cross sections are too small to be amenable to 
measurements. Thus to get the larger cross section (i.e. counting rates) forward 
angles experiments (where q is smaller) with beam energies around 1 GeV appear 
to be favourable. 

In Fig. [5] we show the dependence of the total cross section on piab for the 
reactions studied in Figs. [H] and [T] We note that in both the cases, as beam 
energy increases beyond the threshold, the cross sections first increases rapidly 
and after peaking around a given value starts decreasing slowly. The beam 
momentum where the cross section peaks is around 1.0 GeV. This observation 
was also made in the DWIA study of this reaction in Ref. [1^. This further 
highlights the point made in Fig. [SI 

Next we examine the role of the bound state relativistic effects on the cross 
sections. We would like to point out that since our calculations have been 
performed in momentum space, there is no need to introduce any local approxi- 
mation to the propagators appearing in the production amplitudes — most of the 
non-relativistic calculations necessarily make such approximation. In Fig. IIOI 
we show the role of the lower component of the Dirac bound state wave func- 
tions (g) on the angular distributions of the (7r+,if+) reaction on the ^^C and 
targets for the transitions as depicted in this figure. We see that effect of 
g is more pronounced at larger angles. It should, however, be noted that the 
results with upper components only can not be directly equated with those of 
the conventional approaches that employ the Schodinger equation to describe 
the bound state wave functions. One requires the matrix elements and partic- 
ularly the operator that is given in terms of the Dirac matrices to undergo a 
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Figure 9: The total cross section for the [(Pg^j)^ > (''i/2)a] transition in the tt'^+^'^C 
+^\C reaction (left panel), and the [(c(^y2)jVi (^i/2)a] transition in the 7r++*''Ca 
+^8Ca reaction (right panel) as a function of the incident pion momentum piab- 
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Figure 10: The differential cross section for the same reaction as in the left panel of Fig. |6]and 
for the [(ff^^j)^' {/7/2)a] configuration in the 7r++**^Y — > +^^Y reaction. Shown are the 
results including the full relativistic spinors (solid lines) and without the lower components 
(dotted lines) as a function of the kaon angle 6. 



non-relativistic reduction. Such a reduction has been carried out in Refs 55l.l5q 
where it is shown that for the hypernuclear production by the (7, K^) reaction, 
the difference between the relativistic and non-relativistic cross sections could 
be about 20% even in the forward directions. 



5.1.1. Spectral distributions for ^\ V and ^% Y hypernuclei 

In Figs. [TT] and [T21 we show the comparison of the calculated spectral dis- 
tributions for the ^AV and hypernuclei produced in the (7r+, i4r+) reaction 
on ^^V and ^^Y targets at the beam momentum of 1.05 GeV. In drawing the 
smooth spectral distributions, each level is a convolution with the appropriate 
Gaussian width F depending on its character in the experimental data. We used 
only a single Gaussian for all the levels in contrast to Ref . [3l , where more than 
one Gaussian were used for some levels. 

In ^IV case the series of levels are obtained by the following configura- 
tions and binding energies (i^bind): [{fyl)N, (si/2)a] (3", i^bind = 19.75 MeV), 
[{fr/\)N,PA] (4+, i^bind = 11.75 MeV), and [(/f/Jjv, ^a] (S", £;bind = 3.75 MeV). 
These states make the largest contributions to the corresponding cross sec- 
tions. On the other hand, the levels in are obtained with the config- 
uration where the neutron hole state corresponds to the ^9/2 orbit with the 
A in the Os, Op, Od and 0/ orbitals. The configurations of these levels are: 
[(5-/2)Ar,(0s)A](4+,Sbind = 23.6MeV), [(5-/2)^, (0p)a](5-, i;bind = 16.5 MeV), 
[{g-/\)N, (0d)A](6+,£;bind = lO.OMeV), [{g-/\)N, (0/)a](7-, £;bind = 2.3 MeV). 

We see that for both nuclei the calculated spectral distributions reproduce 
the overall global trends of the data of Ref. [l3| reasonably well. Our calculations 
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reproduce the experimentally observed three and four clear peaks in the ^IV and 
spectra, respectively. However, looking more closely we note that the theory 
overestimates the experimental cross sections for the (Os)a orbitals somewhat 
in both cases. For these states, the corresponding binding energies (i?bind) are 
quite large and for these strongly mismatched cases, distortion effects could play 
a more significant role. 

In both cases, one observed some minor peaks which fill up the gaps between 
the major peaks. These peaks correspond to core excited states and the states 
corresponding to other, less dominant, members of the configurations mentioned 
above. We have not made any attempt to fit to these peaks in this exploratory 
study. Clearly, a quantitative description of the spectral shapes of the heavier 
mass A hypernuclei requires proper consideration of the core excitation and 
mixing of states of different parity. 

5.2. Effects of initial and final state interactions 

The calculations presented thus far have been done in a plane wave ap- 
proximation, where the pion-nucleus and kaon-nucleus interactions are ignored. 
While the essential features of the high momentum transfer reaction (7r+,i^+) 
can be understood in this approach, the nuclear interactions may have some 
consequences. They produce both absorptive and dispersive effects Ref. (20| . 
However, for the large incident energies considered in this calculation, the ab- 
sorption effects are likely to be the more important. 

One should also note that in the distorted wave treatment the continuum 
wave functions are no longer associated with sharp momenta but are states with 
a momentum distribution, see Eqs. (|12p . and This leads to a redistribution 
of the momentum transfer differently from what is allowed in the plane wave 
approximation. It could shift the sensitivity of the model to even lower momenta 
leading to enhanced cross sections. The competition between this effect and 
the absorption effect would ultimately decide the role of distortions in these 
reactions. 

We have estimated the role of the distortion effects in the initial channel on 
the absolute magnitude of the cross sections. For this purpose, we make use of 
the eikonal approximation as discussed in section Pf. II Cross sections have been 
calculated in both PW and eikonal approximation at one most forward angle for 
i2C(7r+, '^°Ca(7r+, if+)fCa and ^^Y{Tr+ , K+)fY reactions. The results 

are shown in table [SI 

We see that the distortion effects in the incident channel lead to the reduction 
of the peak cross sections by factors of 1.2 to 12 as the target mass varies from 12 
to 89. Our result for the '^'^C target is in contrast to that of the non-relativistic 
impulse model of Ref. , where the distortion effects were found to reduce the 
differential cross section by about an order of magnitude for this target. The 
distortion effects play an increasingly important role with increasing mass of the 
target nucleus. This result is in agreement with those of Refs 23], and (Hij . 
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Table 6: The differential cross sections at selected angles for the plane wave approximation 
and the distorted wave calculations for the pion in the eikonal approximation. 



nucleus 


transition 


da/dn [nb/sr] 
plane wave eikonal 


i^C [Ok = 1°) 


np3/2 Asi/2 
npz/2 ^P3/2 


16.55 13.21 
10.16 7.10 


40Ca [Ok = 4°) 


nd3/2 Asi/2 


15.55 3.18 


89 Y {Ok = 4°) 


7109/2 Asi/2 


7.43 0.642 



6. Summary and conclusions 

In summary, wc studied the A{tt^ , K~^)aA reaction on ^^C, ^°Ca, '''^V and 
targets within a fully covariant model, where in the initial collision of inci- 
dent pion with one of the target nucleons the A/'* (1650), A/'* (1710), and A/'* (1720) 
intermediate baryon resonance states are excited which subsequently propagate 
and decay into the relevant channel. We have retained the full field theoretic 
forms of various interaction vertices and obtained the baryon bound states by 
solving the Dirac equation with appropriate scalar and vector potentials. The 
vertex constants were taken to be the same as those determined in previous 
studies. We have ignored, for the time being, the distortion effects in the inci- 
dent and outgoing channels as our main aim in this paper has been to establish 
a fully covariant model for this reaction which has so far been described only 
within a non-relativistic distorted wave impulse approximation picture. The 
plane wave approximation facilitates the application of this novel approach to 
very many cases without requiring lengthy and cumbersome computations which 
are necessarily involved in the distorted wave methods. 

We find that excitations of iV*(1650), and TV* (1710) resonant states domi- 
nate the cross sections for the (tt^, if^) reaction for beam energies below 2 GeV. 
Our model describes well the shapes of the experimental angular distributions 
for two states which corresponds to the prominent peaks in the -'^^C spectrum. 
The total cross sections show a substantial dependence on the beam energy with 
a distinct peak around 1 GeV for reactions on both ^^C and ^°Ca targets. The 
differential cross sections peak near zero degrees for both light mass as well as 
heavier targets. Thus, measurements at forward angles and at beam energies 
around 1.0 GeV are expected to have high yields. 

The characteristic bump structures reflecting the A major shell orbits as 
seen in the spectra of the "''aV and hypernuclei were reproduced reasonably 
well by our model. However, for a quantitative description of the spectra, more 
realistic calculations with configuration mixed A particle, neutron hole states 
are required. Mixing of the different parity states within a given particle-hole 
configuration may also be needed. 

Our work shows that it is indeed feasible to have a fully covariant description 
of the hypernuclear production via the (tt^, K~^) reaction. In future studies, the 
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structure aspects of the model should be improved. Also, alternative covariant 



bound state wave functions [1^, |57| should be used in order to check the validity 
of the method used in this work to determine them. Distortion effects in the 
initial and final channels are required in order to have more reliable predictions 
of the absolute magnitudes of the cross sections. 
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